Mon. Not. R. Astron. Soc. 000, 1 (2006) Printed 4 February 2008 (MN KTpX style file v2.2) 



Current models of the observable consequences of cosmic 
reionization and their detectability 

Ilian T. Iliev 1 *, Garrelt Mellema 2 , Ue-Li Pen 1 , J. Richard Bond 1 , 
Paul R. Shapiro" 



O 
O 

(N 
-i— > 
O 

O 

m 
(N 

CM 
> 

ON 

o 

(N 

o 
o 

6 



1 Canadian Institute for Theoretical Astrophysics, University of Toronto, 60 St. George Street, Toronto, ON M5S 3H8, Canada 

2 Stockholm Observatory, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden 

3 Department of Astronomy, University of Texas, Austin, TX 78712-1083, U.S.A. 



4 February 2008 



ABSTRACT 

A number of large current experiments aim to detect the signatures of the Cosmic Reion- 
ization at redshifts z > 6. Their success depends crucially on understanding the character 
of the reionization process and its observable consequences and designing the best strate- 
gies to use. We use large-scale simulations of cosmic reionization to evaluate the reionization 
signatures at redshifted 21 -cm and small-scale CMB anisotropies in the best current model 
for the background universe, with fundamental cosmological parameters given by WMAP 3- 
year results (WMAP3). We find that the optimal frequency range for observing the "global 
step" of the 21-cm emission is 120-150 MHz, while statistical studies should aim at 140- 
160 MHz, observable by GMRT. Some strongly-nongaussian brightness features should be 
detectable at frequencies up to ~ 190 MHz. In terms of sensitivity-signal trade-off relatively 
low resolutions, corresponding to beams of at least a few arcminutes, are preferable. The CMB 
anisotropy signal from the kinetic Sunyaev-Zel'dovich effect from reionized patches peaks at 
tens of fiK at arcminute scales and has an rms of ~ lfiK, and should be observable by the 
Atacama Cosmology Telescope and the South Pole Telescope. We discuss the various obser- 
vational issues and the uncertainties involved, mostly related to the poorly-known reionization 
parameters and, to a lesser extend, to the uncertainties in the background cosmology. 

Key words: H II regions — high-redshift — galaxies: formation — intergalactic medium — 
cosmology: theory — radiative transfer — methods: numerical 



1 INTRODUCTION 

The Epoch of Reionization and the preceding Cosmic Dark Ages, 
from recombination at z ~ 1100 to z ~ 6 include the for- 
mation of the first nonlinear cosmological structures, first stars, 
QSO's and the emergence of the Cosmic Web as we know it to- 
day. Yet it still remains poorly understood. This is mostly due to 
the scarcity of direct observations, resulting in weak constraints 
on the theoretical models. However, this situation is set to im- 
prove markedly in the coming years, with the construction of a 
number of new observational facilities, particularly for detection of 
the redshifted 21-cm line of hydrogerQ and the kinetic Sunyaev- 
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Zel'dovich (kSZ) effecfl. Simulations of that epoch, required in 
order to make reliable predictions for such observations, are diffi- 
cult and computationally-intensive. Recently we presented a set of 
large-scale, high -resolution radi a tive transfer simulations of cosmic 
reionization ( Iliev et al.l [20063 : iMellema eTall [2006b : Iliev et all 
l2007ah . These simulations were the first ones which were suf- 
ficiently large to reliably capture the characteristic scales of the 
reionization process. This allowed us to derive the first realistic 
predictions for the reionization observables, in particular the differ- 
ent si gnatures in the redsh ifted 21-cm emission of neutral hydro- 
gen dMellema et a l. 2006b) and the imprint of the ionized patches 
on small-scale CMB temperature anisotropies through the k inetic 
Sunyaev-Zel'dovich (kSZ) effect dlliev et al J l2006cl . 1 2007bh . We 
also proposed an approach to use the obtained 21-cm maps to de- 
rive the Thomson optical depth fluctuations due to reionization 
(Holde ret~aLll2006l) and derived the CMB polarization signatures 
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of patchy reionization jPore et ai]|2007b . We studied the effects of 
varying ionizing source efficiencies and sub-grid gas cl umping on 
these observed signals. All of these calculations except iDore et al.l 
d2007h used a particular set of cosmological parameters, based 
on the best-fit first-year WMAP results, hereafter WMAP1, with 
the following cosmological para meters: (fl m , ^a, flfr, h, erg, n) = 
(0.27,0.73,0.044,0.7,0.9,1) dSpergel & et aT]|2003l) . Here fi m , 
Qa, and Qt are the total matter, vacuum, and baryonic densities 
in units of the critical density, erg is the rms density fluctuations 
extrapolated to the present on the scale of 8/i _1 Mpc according to 
the linear perturbation theory, and n is the index of the primordial 
power spectrum of density fluctuations. 

Recently t he 3-y ear WMAP results were published 
dSpergel & etall |2006h , hereafter WMAP3, which pre- 
sented an updated, and fairly different best-fit cosmology: 
(tl m ,tl A ,tl b ,h,a s ,n) = (0.24,0.76,0.042,0.73,0.74,0.95). 
We also note that other current measurements of these parameters 
based on e.g. other CMB experiments, supernovae, large-scale 
structure, clusters, and Ly-a forest tend to give slightly different 
valu es, either on their own or in combination with the WMAP dat a 
(e.g. ISpergel & et al]|2006l : ISeliak et al.ll20ol : lYao & et al.ll2006h . 
In particular, they tend to yield higher value of erg than WMAP3 
alone (but still well below the WMAP1 value), at erg ~ 0.8 - 0.85. 
For example, recent results using all of the CMB data derives 
as = 0.79, whil e combining with the large-scale structure data 
yields a a = 0.81 l lKuo & et alfeOOrj) . 

In this work we evaluate the detectability of reionization at 
radio wavelength observations of redshifted 21 -cm line of hydro- 
gen and CMB anisotropies from kSZ. We derive a variety of ob- 
servational signatures and discuss their detectability with a num- 
ber of current and near-future experiments and discuss the related 
uncertainties due to various poorly-known parameters. We give 
special attention to matching the parameters of the observations 
(beamsizes, bandwidths, frequencies, observational strategies) to 
the character of the reionization features in order to maximize their 
detectability. We also compare our predicted signals to the expected 
sensitivities for several current experiments. 

For the benefit of the reader, whenever possible we compare 
our current results with our previous predictions done in the frame- 
work of the WMAP1 cosmology. The major difference between 
the WMAP1 and WMAP3 cosmologies is the overall amplitude of 
the power spectrum, expressed here in terms of erg, but the mod- 
els also have slightly different spectral shapes, with the low erg one 
having a red tilt, n „ — 1 = —0.05. We have previously shown 
dAlvarez et ai]l2006l :l Iliev et al. 2007a) that these changes result in 
structure formation being delayed in the WMAP3 universe relative 
to the WMAP1 universe, so the epoch of reionization is shifted to 
lower redshifts. In particular, if source halos of a given mass are as- 
sumed to have released ionizing photons with the same efficiency 
in either case, then reionization for WMAP3 is predicted to have 
occurred at (1 + z)-values which are roughly 1.3-1.4 times smaller 
than for WMAP1. The predicted electron-scattering optical depth 
of the IGM accumulated since the beginning of the EOR would 
have then been smaller for WMAP3 than for WMAP1 by a factor 
of (1.3— 1.4) 3//2 ~ 1.5— 1.7, just as the observations of large-angle 
fluctuations in the CMB polarization require. This means that the 
ionizing efficiency per collapsed baryon required to make reioniza- 
tion early enough to explain the value of r cs reported for WMAP1 
and WMAP3 are nearly the same. 

This delay of reionization can be understood in terms of the 
density fluctuations at the scales relevant to reionization as follows. 
Let us denote the rms linear amplitudes on the top hat smooth- 



ing scales of 0.1 h Mpc and 0.01 h 1 Mpc by cro.i and co.oi, re- 
spectively. The top hat scale 0.1 /i _1 Mpc corresponds to a mass 
2.7 x 10 8 /i _1 M for the low erg case, and 3.1 x 10 8 /i _1 M Q for 
the high erg case, the slight difference being due to the differing Q m 
and h. Of course 0.01 h~~ 1 Mpc corresponds to masses 3 orders of 
magnitude smaller, 2.7 x 10 5 /i _1 M© and 3.1 x 10 5 /i _1 M , for 
WMAP3 and WMAP1. Thus, cro.i and oo.oi roughly correspond to 
the scales of the dwarf galaxies and minihaloes, respectively. The 
shape difference in the low and high erg cases is encoded in the 
ratios aui/og, 6.1 and 6.6, and <7o.oi/<7g, 10.1 and 11.0, that is, 
not negligible but not that large relative to the 20% decrease in erg. 
A reasonable indication of when structure on scale R formed at 
high redshift is 1 + z R w 1.3cr fl n~ ' 23 , where [a/D(z)] « O^ 23 
for jz 3> 1. Here D(z) is the linear growth factor from red- 
shift z to the present. For the minihalo scale R = 0.01 /i _1 Mpc, 
zo.oi ~ 12.5 and 16.4, respectively. For the dwarf scale 20.1 ~ 
7.1 and 9.4, respectively, in reasonable accord with the computed 
overlap redshifts f rom our inhomogeneous reionization simulations 
dlliev et alj|2007ij, and Table Q] below); the 50% reionization red- 
shifts bracket the minihalo and dwarf structure formation redshifts. 
The uniform reionization Thompson depth r to a reionization red- 
shift z re i is t = 0.085 [( 1 + z re i ) / 1 1] 3 ' 2 for the low erg model and 
t = 0.080[(l+2r e i)/ll] 3/2 for the high erg model. When the z 50% 
values are substituted, there is rough agreement with the r cs in Ta- 
ble[TJ The scaling of the Thompson depth r would be (1 + 2_r) 3 ^ 2 , 
about 1.4, roughly co nsistent with the sa 1.3 ratio we determine , 
and with the results in lAlvarez et al] {2006) an d Iliev etal. (2007a). 

Some of the observable implications of this delay in the for- 
mation of structures are fairly straightforward. For example, the 
decrease in the mean spatially-averaged redshift ed 21 -cm signal as 
the IGM reionizes, referred to as "global step" l lShaver et al.ll 19991) 
will occur at lower redshift, higher frequency, in WMAP3 cosmol- 
ogy. It should thus become a bit easier to observe than previously 
thought, due to the lower foregrounds and higher sensitivity at high 
frequencies. However, other consequences of the new cosmology 
framework are less obvious and have to be evaluated with care. 

This paper is organized as follows. In §[2] we briefly describe 
our simulations. In § [3] we present our predictions for the red- 
shifted 21 -cm signals and discuss their observability with current 
and planned radio arrays. In §[4] we evaluate the patchy kSZ signal 
and its observability with ACT and SPT telescopes. Our conclu- 
sions are summarized in §[5] 



2 SIMULATIONS 

Our simulations were performed using a combination of two very 
efficient comput ational tools, a c osmological particle-mesh code 
called PMFAST dMerz et alj|2005h for following the structure for- 
mation, whose outputs are then post-processed using our radia- 
tive transfer and non- equilibrium chemistry code called C 2 -Ray 
dMellema etai] |2006a). Our simulations, p arameters and method- 
ology were discussed in llliev et al. 12006b): Mellema et alJ d2~006b) 
and llliev et al. Detailed tests of our radiative transfer 

method were presented in Mell ema et al. | d2006al) and llliev et all 
(2006a). The simulations considered in this work are summarized 
in Table [T] along with the basic characteristics of their reioniza- 
tion historiefl The parameter / 7 characterizes the emissivity of 
the ionizing sources - how many ionizing photons per gas atom in 

3 WMAP1 cases listed here were first presented in <Mellema et alj2006bl) . 
with predictions of 21 -cm background and kSZ effect from patchy reion- 
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Table 1. Simulation parameters and global reionization history results for 
the simulations used in this work. 



?S - Tomb 2>\iA w T*nHi{z) 





f250 


f'250C 


mesh 


203 3 


203 3 


box size [ h _1 Mpc] 


100 


100 


ft 


250 


250 


'-'subgrid 


1 


C(z) 


z 50% 


8.9 


8.3 


-^overlap 


7.5 


6.6 


7~cs 


0.082 


0.076 



the (resolved) halos are produced and manage to escape from the 
host halo within ~ 20 Myr, which is the time between two con- 
secutive density slices, equal to two radiative transfer timesteps. 
Both of our WMAP3 simulations utilize / 7 = 250, as this value 
yields final H II region overlap at z ov ~ 6.5 — 7.5 in agreement 
with the current observational constraints. The corresponding in- 
tegrated Thomson scattering optical depths, r cs are also in agree- 
ment, within 1-ct of the WMAP3 derived value, r cs = 0.09 ± 0.03, 
although they are a bit lower than the central value. Much lower 
ionizing efficiencies for the sources than these will result in too 
late an overlap, violating the available observational constraints, 
while much higher ones would result in an early-reionization sce- 
nario, again possibly in conflict with observations. The correspond- 
ing WMAP1 simulations with the same ionizing source efficiencies 
of / 7 = 250 yielded r cs = 0.098 — 0.130, outside of the nominal 
1-ct WMAP1 -derived range r eB = 0.17 ± 0.04. As we have previ- 
ously shown dlliev et al.f2 007a), the presence of low-mass ionizing 
sources [abs ent here, but resol ved in the smaller-box simulations 
presented in llliev et al. (2007a)] increases the total optical depth, 
and can easily bring it into agreement with any value within the 
WMAP3 (or WMAP1, respectively) 1-ct range. In the same previ- 
ous work we showed that despite this potentially dramatic effect on 
the integrated Thomson optical depth, the presence of small sources 
has only modest effects on the large-scale geometry of reioniza- 
tion, because most of these smaller sources were strongly clustered 
and as a consequence become strongly suppressed during the later 
stages of reionization through Jeans-mass filtering in the ionized 
regions. 



3 21-CM EMISSION 

The differential brightness temperature with respect to the CMB 
of the redshifted 21 -cm emission is determined by the density of 
neutral hydrogen, pui, and its spin temperature, Ts. It is given by 

Ts — Tcmb 



5T, 



1 + z 



-(1 



1 + z 
= 28.5 mK 



32nT s H(z) 



(i) 

'0.24\ 1 /2 



/0.24V 
V Mm ) 



,0.042 0.73, 

dFieldlll959h . where z is the redshift, Tcmb is the temperature of 
the CMB radiation at that redshift, r is the corresponding 21-cm 
optical depth, Ao = 21.16 cm is the rest-frame wavelength of 
the line, A\q = 2.85 x 10 _15 s _1 is the Einstein A-coefficient, 
T» = 0.068 K corresponds to the energy difference between the 
two levels, 1 + S — nm/{riH) is the mean number density of neu- 
tral hydrogen in units of the mean number density of hydrogen at 
redshift z, 

{n H ){z) = —1 + z 

HHm p 

= 1.909 xl0- 7 cm- 3 (^) (1 + z) 3 , (2) 

with fijj = 1.32 the corresponding mean molecular weight (as- 
suming 32% He abundance by mass), and H{z) is the redshift- 
dependent Hubble constant, 

H(z) = # [fim(l + zf + n k (l + Z f + fl A ] 1/2 



= H E(z) ^H nK 2 (l + z 



\3/2 



(3) 



ization from those cases presented in Mellema et al. ( 2006b) and llliev et alj 
l2007bl) . respectively. 



where Ho is its value at present, and the last approximation is valid 
for z S> 1. Throughout this work we assume that all of the neutral 
IGM gas is Ly-a-pumped and sufficiently hot (due to e.g. a small 
amount of X-ray heating) above the CMB temperature and is thus 
seen in emission. These assumptions are g enerally well-justified , 
except possibly at the earliest times (see e.g. lFurlanetto et al.l200rj . 
and references therein). 



3.1 The progress of reionization: global view 

The images in Figure Q] show the progress of reionization as seen 
at 21-cm emission vs. the observed frequency for run f250. The 
technique we used to produce them was described in detail in 
iMellema et all {2006b) and is similar to the standard method for 
making light-cone images. In short, it involves continuous in- 
terpolation between the single-redshift numerical outputs in fre- 
quency/redshift space. The slices are cut at an oblique angle so as 
to avoid repetition of the same structures along the same line-of- 
sight. Redshift-space distortions due to peculiar bulk velocities are 
included. The top image shows the decimal logarithm of the differ- 
ential brightness temperature at the full resolution of our simulation 
data, approximately 0.25 arcmin in angle and 30 kHz in frequency. 
The ionization starts from the highest density peaks, which is where 
the very first galaxies form. These high peaks are strongly clustered 
at high redshifts, which results in a quick local percolation of the 
ionized regions. As a result, by z ~ 9 (fobs ~ 140 MHz) already a 
few fairly large ionized regions form, each of size ~ 10 Mpc. These 
continue to grow and overlap until there is only one topologically- 
connected H II region in our computational volume, which for this 
particular simulation occurs at z ~ 7 (i^obs ~ 180 MHz). However, 
even at this time quite large, tens of Mpc across, neutral regions still 
remain. They are gradually ionized as time goes on, but some of 
them persist until the very end of our simulation. The uncertaities 
in the reionization parameters (e.g. source efficiencies, gas clump- 
ing) result in moderate variations in the typical sizes of the ionized 
and neutral regions, but do not change the basic characteristics of 
the reionization process. 

The bottom image in Figure[T]shows the same data as the top, 



4 I. T. Iliev, et al. 




50- 

.a 40 

3 20- 
£ 10 ■ 



J . 1 I. 




2,26 
3,69 



180 



■■■'»■»■ 



170 



150 



140 



160 
V (MHz) 

i 1 ' J J J ' ■ ■ ■ [ ■ ■ ■ L ' / 1 \,l * J L J. ' J , ' ' ' ,' 1 ,'■ ' 



-n- 



i i | i i i 

180 



170 



M I I I I I 



-ri- 



160 
V (MHz.) 



1 \ T 

150 



t i i i r 



122.7 
12.5 

2.35 

-7.83 
1-18.0 



140 



Figure 1. Position-redshift slices from the f250C simulation. These slices illustrate the large-scale geometry of reionization and the significant local variations 
in reionization history as seen at redshifted 21 -cm line. Observationally they correspond to slices through an image-frequency volume. The top image shows 
the decimal log of the differential brightness temperature at the full grid resolution, while the bottom image shows the same T;, data, but in linear scale and 
smoothed with a compensated Gaussian beam of 3' and (tophat) bandwidth of 0.2 MHz, roughly corresponding to the expected parameters for LOFAR. The 
spatial scale is given in comoving units (top) or (approximate) angle on the sky in arcminutes (bottom). 




^ obs [MHz] 

Figure 2. Evolution of the mean mass-weighted ionized fraction, x m , (top), 
mean flux in fi Jy (assuming a circular beam with an angular diameter of 
10') (middle) and the mean differential brightness temperature, STf,, in mK 
vs. observed frequency for simulations f250 (solid, blue), and f250C (short- 
dashed, red). 



but as would be seen by a radio interferometer array assuming per- 
fect foreground removal. To obtain it we smoothed the data with 
a compensated Gaussian beam with FWHM of 3' and bandwidth 
of 0.2 MHz, similar to the expected parameters for the LOFAR ar- 
ray. Unless otherwise stated, throughout this paper we use a com- 
pensated Gaussian beam. We recognise that a compensated Gaus- 



sian is not perfect match to the actual interferometer beam, how- 
ever it captures its essential properties. In particular, it has zero 
mea n, and negative troug hs at the side of the central peak (see 
e.g. iMellema et aUl2006bh . As a result, ionized regions show as 
negative differential brightness temperature regions if they are sur- 
rounded by neutral volumes. If an ionized (or a neutral) region is 
much larger than the smoothing scale, the resulting signal would be 
close to zero. A direct comparison between the two images shows 
that all the main structures clearly remain also in the smoothed 
image, indicating that LOFAR and other similar interferometers 
would have sufficient resolution to determine the large-scale reion- 
ization morphology to a reasonable accuracy throughout most of 
the reionization history. However, at the earliest stages of reioniza- 
tion some of the existing H II regions are barely seen, or not at all, 
since the beam smoothing either merges them with other nearby 
ionized regions, or simply smoothes them away. This is a conse- 
quence of the small sizes of these early H II regions, which puts the 
majority of them below the beam resolution. At frequencies higher 
than ~ 140 MHz the ionized bubbles become large enough to be 
above the smoothing scale and thus all major structures become 
visible. Higher ionizing source efficiencies and/or lower gas clump- 
ing would yield somewhat larger H II regions, thus they would be- 
come visible slightly earlier. We also note that the Cosmic Dark 
Ages and the early stages of reionization might be still observable 
through other sources of 21 -cm fluctuations which we do not con- 
sider here. These include e.g. the 21-cm emis s ion and absorption 
by cosmological mini haloes Jlliev et alj[2003 : Furlanctto & Loeb 
20021 : llliev etldll 20031), by shock-heated IGM iFurlanetto & Loeb 
20041 : IShapiro et alj|2006|). or due to inhomogeneous early ba ck- 
ground s in Lv-q |Barkana & LoeblfeOQSl : IChuzhov et alj|2006l) or 
X-rays jPritchard & Furlanettoll200q) . We note that there are bright 
and potentially observable features even for v a i, B > 180 MHz. 
This is important for the planned observations since at such high 
frequencies both the foregrounds and the Radio Frequency Inter- 
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ference (RFI) are substantially lower, while at the same time the 
array sensitivities improve. We discuss these issues in more detail 
in §[331 

In Figure [2] we show the evolution of the mean (i.e. aver- 
aged over a large volume) mass-weighted ionization fraction, x m , 
and the mean differential brightness temperature as seen at the 
observer, 5Tb(v bs)- As the universe steadily becomes ever more 
ionized, the mean 21 -cm differential brightness temperature natu- 
rally decreases. However, the signal persists at a non-trivial level 
(a few mK or more) until quite late,up to frequencies of 150-170 
MHz. Similarly to our previous results which used the WMAP1 pa- 
rameters, the_^gj^bal ster£_from mostly-neutral to mostly-ionized 
medium jShaveretal.1ll999h turns out to be rather gradual, with 
the signal decreasing by ~ 20 mK over 20-30 MHz (and some- 
what more gradual for case f250C than for f250). The differen- 
tial brightness temperature scales with the reionization redshift as 
5T b oc (1 + z) l/2 . Thus, a delay of reionization from WMAP1 to 
WMAP3 model by 1.3 in redshift corresponds to an expected 8Tb 
decrease by a factor of 1.14, roughly as found in our simulations. 

Analytical models of th e globally-averaged 21 -cm signal 
dSethil 120051 : rFurlanettdl2o61. e.g .) predict evolutions which are 
in rough ag reement with ones we find. The models considered in 
ISethil (2005) yielded a somewhat faster evolution of the global sig- 
nal, while the wider range of models presented in lFurlanettd ( 2006) 
included cases of both faster and of more gradual evolution, de- 
pending on the details of the evolution of the Ly-a and X-ray back- 
grounds. 

3.2 The statistical measures 

An alternative approach to take would be a statistical one, through 
detection of the fluctuations of the emis sion around its mean 
value (see Scott & Rees 1990; Madau et al. 



1997l:lllievetal.ll2002l; 



Morales & Hewittl |2004 iMellema etal] 



Zaldarriaga et alj 120041; 
2006b; iFurlanetto et al.ll2006l and references therein). These fluc- 
tuations are due to a combination of the reionization patchiness 
and the variations of the underlying density field. In Figure [3] we 
show the rms of the 21-cm emission fluctuations derived from our 
simulations. The top panel shows the evolution of the differential 
brightness temperature rms, 5Tb, rm s, for three choices for the beam 
size (in arcmin; using compensated Gaussian beam) and the band- 
width (in MHz), from top to bottom, (A#beam, Atw) = (0.1, 1) 
(roughly as expected for the SKA compact core), (0.2, 3) (LOFAR) 
and (0.4, 6) (GMRT, MWA). The middle panel shows the evolution 
of the corresponding fluxes, given by 



8F{v ohs ) = 



-kB5Tb(v bs)Afl, 



(4) 



where AS1 = ir(8/2) 2 is the 3D angle subtended by the beam 
for a circular beam with FWHM of #beam- On the bottom panel we 
show the 21-cm signal as a fraction of the dominant foreground, the 
Galactic synchrotron emission. For the larger beams/bandwidths 
the 5T bl rms peak is at z ~ 8.8 (v ohs ~ 145 MHz) for f250 and at 
z ~ 8.3 {v ohs ~ 152 - 153 MHz to for f250C, close to the time 
at which 50% of the mass is ionized, a s was the case also for our 
WMAP1 results [Mel lema et afl f2006b). However, for the higher 
resolution, corresponding to a 1' beam, the peak of the fluctuations 
moves to noticeably earlier times/lower frequencies, to ~ 137 MHz 
(x m = 0.26) for f250 and to ~ 143 MHz (x m = 0.26) for f250C. 
This high-resolution case differs from the rest because the scales 
probed by such a small beam/bandwidth combination are generally 
below the characteristic bubble size. The exception is at early times, 
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Figure 3. (a)(top) rms fluctuations of the differential brightness tempera- 
ture, <57b rms = (ST 2 ) 1 / 2 vs. observed frequency, v b s for f250 (blue) 
and f'250C (red) for beam sizes and bandwidths ASbcam, Afb w ) =(1,0.1 
MHz) (dotted), (3', 0.2 MHz) (solid) and (6', 0.4 MHz) (dashed); 
(b)(middle) Fluxes in fijy corresponding to the differential brightness tem- 
perature fluctuations in (a), same notation as in (a); and (c)(bottom) Ratio 
of <5T b rms to the mean Galactic synchrotron foreground, assumed to be 



<57b, ga i = 300K(i//150MHz)- 



, same notation as in (a). 



when ionized bubbles are still small on average, and thus match 
better the smaller beam-size, which is reflected in the fluctuation 
peak moving to earlier times. 

We note that the flux fluctuations peak somewhat later than the 
corresponding temperature fluctuations. The position of the peak 
flux in redshift/frequency space is largely independent of the res- 
olution employed, and is at ~ 148 — 150 MHz for f250 and at 
~ 155 — 157 MHz for f250C. As the beam-size and bandwidth 
increase from 1' to 3' and then to 6', the temperature fluctuations 
decrease, albeit only by a modest amount. For example, the peak 
fluctuations for (A#beam, Aia,w) =(6',0.4 MHz) are only a factor 
of ~ 2 lower (4 mK vs. 8 mK), while the corresponding flux in- 
creases by over an order of magnitude, indicating that it would be 
optimal to observe at relatively large scales, where we maximize 
the sensitivity without sacrificing much of the signal. 

Higher source efficiency and lower gas clumping would re- 
sult in the bubbles at the same redshift being a bit larger in typ- 
ical size. This would move the peak position to a bit lower fre- 
quency. It would also make it slightly higher for the larger beam 
sizes since larger bubble sizes at the half-ionized point would match 
these beamsizes better. 

The 21-cm temperature fluctuations as fraction of the domi- 
nant foreground, the Galactic synchrotron emission (bottom panel) 
peak even later than the flux fluctuations. This is due to the broad 
peak of the 21-cm emission and the steep decline of this foreground 
at higher frequencies, which combine to push the peak to later 
times/higher frequencies, at v ^ s up to 160-165 MHz. The signal 
is strongly dominated by the foregrounds at all times, but up to an 
order of magnitude could be gained for observations aimed close to 



6 I. T. Iliev, et al. 



0.6 
0.4 

0.2 

0.1 

0.6 
0.4 

0.2 

0.1 



_ f250 

z=12 


z = 8.8 


ii i i nun i i mi 

z = 8.5 


1 1 1 llllll 1 1 llllll \ 
z = 7.5 



0.1 



1 



10 0.1 



1 



10 



< 
< 



0.6 
0.4 

0.2 

0.1 

0.6 
0.4 

0.2 

0.1 



_ f250C 

z=12 


z = 8.3 


llllll 1 1 llllll 1 1 llllll ; 

z = 7.9 


llllll 1 1 llllll 1 1 llllll \ 
z = 6.6 



0.1 



1 



10 0.1 



1 



10 



k [h/Mpc] k [h/Mpc] 

Figure 5. The variance, A of the 3D power spectra of the neutral gas density A p ^hi (dotted), and the total density A p (solid), normalized to the total, for 
simulations f250 (left) and f250C (right). The redshifts are chosen as follows (in decreasing numerical order): early times (z = 12; Tf, = 31.8 mK for f25(), 
% = 33.2 mK for f250C), the redshift at which <5T bjI . ms peaks (z = 8.8, % = 13.1 mK for f250, z = 8.3, % = 11.8 mK for f250C), the redshift at which 
5F IIas peaks (z = 8.5, % = 8.9 mK for f250; z = 7.9, % = 8.2 mK for f250C), and the redshift of overlap (5£ = 0.3 mK for f250; % = 0.2 mK for 
f'250C). 




Figure 4. (a)(top) rms fluctuations of the differential brightness tempera- 
ture, 5T brms = {<5T 2 ) 1//2 vs. beam-size, 5#beam f° r ^50 (solid, blue) 
and f250C (short-dashed, red) at redshifts close to the maximum of the 
fluctuations, as indicated. The bandwidth is changed in proportion to the 
beamsize. (b)(bottom) Flux fluctuations corresponding to the differential 
brightness temperature fluctuations in (a), same notation as in (a). 



the peak ratio, as opposed to earlier or later times. Furthermore, ob- 
serving with an interferometer removes significant part of the fore- 
grounds, due to the differential nature of the measurements. This 
eliminates the uniform component of the foregrounds, leaving only 
its fluctuations, at the level of 1-10% of the total. 



In Figure [4] we plot the differential brightness tempera- 
ture fluctuations and the corresponding fluxes vs. the instrument 
smoothing, as given by the beam size, with the bandwidth changed 
in proportion to the beam size. At small scales (below the typical 
ionized bubble size, a few arcmin or less) the temperature fluctua- 
tions are fairly large (> 6 mK) and dominated by Poisson fluctua- 
tions (e.g. a cell is either ionized or else is neutral). However, the 
corresponding flux is tiny, below 1 /i Jy. Around the typical bubble 
scale (~ 5 — 10') the fluctuations are slightly lower, at ~ 3 — 4 mK, 
but the flux grows strongly, as <5#bcam an d reaches ~ 5fj,Jy at 
SObeam ~ 10'. At even larger angles/bandwidths the ionization 
fluctuations gradually start contributing ever less to the total and 
in the large-scale limit the temperature fluctuations just follow the 
underlying large-scale density fluctuations (multiplied by the mean 
differential brightness temperature). The flux curve should gradu- 
ally flatten and at some point there would be little to be gained by 
a further increase of the observed sky patch since the gain in angle 
is canceled by the decrease of the temperature fluctuations. Thus, 
again there is a trade-off between the signal level and the array sen- 
sitivity. The optimal scale for observations will depend on the best 
sensitivity which could be achieved by a particular radio array. For 
compact arrays the optimal beam size is around 10-20', but could 
be lower than that for sensitive arrays with large collecting area like 
SKA. 

Finally, in Figure [5] we show the three-dimensional power 
spectra of the density and the neutral gas density fluctuations (to 
which the 21 -cm signal is directly proportional) at several illustra- 
tive redshifts. Initially, the neutral density power largely tracks the 
one of the density field, since most of the gas is still neutral. At in- 
termediate and late times the ionization fraction inhomogeneities 
introduce a peak around the characteristic scale of the ionized 
patches (which is k ~ 1 hMpc^ 1 for f250C, and slightly lower 
for f250). Around the time when the fluctuations reach their peak 
the patchiness boosts the power on large scales by factor of ~ 2 
(1.5) for f250 (f250C) compared to the density power spectrum. 



Observability of Cosmic Reionization 7 



x 



1 

0.8 
0.6 

0.4 



~i — i — i — i i i | 1 1 — i — i — i — i i i | 

GMRT 100 hours, zenith 




0.2 



0.1 1 
k [h^Mpc] 

Figure 6. Observability of the 21-cm signal: the 3D power spectrum of the 
neutral hydrogen density, A p jjj, at redshift z = 8.59 (T;, = 16.3 mK) 
with the forecast error bars for 100 hours observation with GMRT vs. 
wavenumber k. We assumed 15 MHz observing bandwidth (the full instan- 
taneous bandwidth of GMRT), T sys = 480 K and T s > T CMB . The 
array configuration is assumed pointed to the zenith, but the sensitivity is 
only weakly dependent on the pointing. 

The signal largely disappears around the time of overlap, since lit- 
tle neutral hydrogen remains, but the power spectra still show the 
characteristic peak, at approximately the same scales. 

Some of the most visible redshifted 21-cm features would be 
the points of maximum departure of the signal from the mean. The 
magnitude of the signal is dependent on the level of smoothing and 
could evolve with redshift. The maxima/minima are roughly inde- 
pendent of redshift, within factor of ~ 2, but there are also some in- 
teresting and nontrivial features. For the case of no smoothing (i.e. 
at full grid resolution, ~ 0.25', 30 kHz) the maximum amplitude 
is quite high, at ~ 100 — 200 mK. Naturally, the beam smoothing 
decreases the amplitude, to ~ 20 — 30 mK for (A#beam, Az^) = 
(1', 0.1 MHz), 10-20 mK for (A0 bcs . m , Av bw ) = (3', 0.2 MHz) 
and to 8-20 mK for (A0 beam , Av bw ) = (6', 0.4 MHz). The intro- 
duction of sub-grid clumping (f250 vs. f250C) results in only minor 
variations here. The absolute value of the minimum is similar to the 
one for the maximum in all cases, but the two still show some dif- 
ferences. For the larger beams/bandwidths the maxima peak around 
the time of 50% ionization by mass, while the absolute values of the 
minima peak noticeably earlier. This is readily understood based on 
the evolution of the typical ionized and neutral region sizes and the 
properties of the compensated Gaussian beam. Early-on the ionized 
regions are small and isolated, surrounded by large neutral patches, 
which yields deep negative minima at the positions of the ionized 
bubbles. At late times the H II regions grow large and the beam- 
smoothed signal inside them is close to zero. The positive maxima, 
on the other hand, are due to the densest neutral regions and are 
sampled by the central maximum of the beam, and thus reach their 
peak later. 

The statistics of these emission peaks is also of considerable 
interest since it shows how common such features are and thus 



what patch of sky one needs to study for a detection. PDFs of 
the 21-cm emission are si milar to the ones w e previously found 
for the WMAP1 cosmology Jlliev et al.l2006bl) . They are consider- 
ably non-Gaussian, especially at late times and smaller smoothing 
scales. In particular, for smoothing scales of ~ 5 — 10 Mpc there 
is an over-abundance of the brightest regions by up to an order of 
magnitude compared to the a Gaussian with the same mean and 
rms. At large scales ( ^ 20 hT x Mpc) and late times the PDFs be- 
come very close to Gaussian. 



3.3 Observability: redshifted 21-cm 

There are several current or upcoming experiments which aim to 
detect the redshifted 21-cm signatures of reionization, including 
LOFAR, MWA, GMRT, PAST/2 1CMA and SKA. Not all details of 
the design and the instruments are yet known, particularly for SKA 
for which even the basic concept is not yet finalized. The only ar- 
rays currently in operation are GMRT and PAST/2 1CMA. Among 
these interferometers, GMRT and LOFAR have the largest collect- 
ing area, at ~ 50, 000 m 2 for GMRT and (effective) ~ 10 5 m 2 for 
LOFAR at 110-200 MHz, and thus in principle they have the best 
sensitivity before the commissioning of SKA. However, the same 
two arrays also have significant interference problems to overcome 
from terrestrial sources of confusion. As an example, in Figure [6] 
we show our predicted 21-cm signal at z = 8.6 (case f250C) along 
with the anticipated GMRT sensitivity to the 3D power spectrum 
for 100 hours of integration. The sensitivity is calculated using 
the baseline distribution as seen from zenith, angle averaged, and 
optimally weighted. Radially, the velocity resolution has been as- 
sumed higher than any scale of interest, which is achievable with 
the recently-developed new software correlator. The error on modes 
are anisotropic for radial versus azimuth and the errors are com- 
bined by weighting each mode by its sum of noise and signal vari- 
ance. We assumed 15 MHz observing bandwidth (the full instanta- 
neous bandwidth of the GMRT correlator, the frequency resolution 
is much better, of order kHz) and T sys — 480 K. The errors calcu- 
lation assumes Gaussian noise and Gaussian source statistics, and a 
7 square degree field of view. The two small-fc bins would probably 
not be observable due to foreground removal. We note that previ- 
ous literature concentrated on two dimensional noise models for 
experiments, and for reference in Figure [7J we show the expected 
band-averaged azimuthal GMRT thermal noise power spectrum for 
bandwidth of 15 MHz, 100 hours of integration, and a single point- 
ing. Noise is estimated using the standard radio telescope noise es- 
timation formulae: 



-Z map — , i 



(5) 



where B is the angular-scale dependent effective beam response 
obtained by summing over all the baselines assuming the tele- 
scope is pointed at the zenith. It changes only slightly when looking 
elsewhere. The beamnoise then is simply the Fourier transform of 
T map . The other two arrays, PAST/21 CM A (~ 10 4 m 2 effective 
area) and MWA (~ 7000 m 2 ) are significantly smaller than either 
LOFAR or GMRT, but are also more compact and would be con- 
structed in areas with very low interference. Thus, they could also 
be quite competitive, at least in terms of detecting 21-cm fluctua- 
tions at very large scales, because of their low spatial resolution. 

Our results indicate that the 21-cm fluctuations in WMAP3 
cosmology peak around u obs — 130 — 170 MHz (compared to 90- 
120 MHz for WMAP1 models), depending on the resolution and 
the detailed reionization parameters (ionizing source efficiencies 
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Figure 7. Observability of the 21-cm signal: forecast angular power spec- 
trum of thermal noise of GMRT for bandwidth of 15 MHz, 100 hours of 
integration, one single pointing. 

and gas clumping). The corresponding 21-cm flux fluctuations, and 
the 21-cm temperature fluctuations as fraction of the foregrounds 
peak at even higher frequencies. These properties significantly fa- 
cilitate the signal detection as compared to the WMAP1 cases, 
since at higher frequencies the detector sensitivities are dramati- 
cally better, while at the same time the foregrounds are lower, by 
roughly factor of ~ (120/170) 2 6 = 2.5. The peak value of the 
differential brightness temperature fluctuations is ~ 4 — 8 mK, de- 
creasing modestly with increasing beam and bandwidth. The cor- 
responding fluxes are much more strongly dependent on the scale 
of observations, ranging from < lfj, Jy for 1' beam and 0.1 MHz 
bandwidth to ~ 8/i Jy for 6' beam and 0.4 MHz bandwidth. This 
clearly argues for relatively large beams, of at least a few arcmin, 
and for maximum collecting area concentrated in the compact core, 
in order to improve the sensitivity while not losing a significant 
fraction of the signal. For example, the expected sensitivity of LO- 
FAR for 1 hour of integration time is ~ 0.6 mjy for bandwidth 
of 1 MHz (and somewhat worse for the virtual core only), while 
its resolution is 3.5"-6" for the whole array and 3'-5' for the core 
(numbers are frequency-dependent). Thus, only the virtual core 
would be useful for EOR observations (although the extended part 
would be important for e.g. better foreground subtraction). Up- 
wards of 100 hours of integration time would be needed for de- 
tection, less if more antennae are placed in the virtual core. 

The importance of aiming at the appropriate frequencies and 
length scales cannot be overstated. As we have shown, around the 
peak of the fluctuations the signal as fraction of the dominant fore- 
ground is up to an order of magnitude larger than it is off the peak. 
Utilizing an appropriate beam, i.e. one which is well-matched to the 
scale of the expected fluctuations and corresponding to sufficiently 
large flux level, given the available sensitivity, is also extremely 
important. Our results indicate that the best beam sizes are of at 
least a few arcminutes, and possibly much larger, of order 5-10' 
or more. Larger ionizing source efficiencies or lower gas clump- 
ing would change the height and position of the peak fluctuations 



slightly, but our main conclusions in terms of the best observation 
strategies and frequency band to observe would remain valid. Once 
more observation data becomes available the reionization scenar- 
ios presented here can be refined to give us better constraints on the 
ionizing sources. 

A nother observationa l strategy, and one of the first to be pro- 
posed ^Shaver et al.l 1 199^) is to aim to detect the "global step" 
over the whole sky which reflects the transition from the fully- 
neut ral universe before reionization to the fully-ionized one after 
(see IShaver et all 1 19991 . for more detailed discussion of the vari- 
ous observational issues). Since this is a global, all-sky signal it 
imposes essentially no requirements in terms of resolution. The 
sensitivity requirements are also relatively modest, since the flux 
corresponding to such a large area on the sky is very large. The 
main difficulty is the foreground subtraction. The global step we 
find is relatively gradual, ~ 20 mK over ~ 20 MHz. This would 
still be readily detectable in absence of foregrounds. How well the 
foregrounds could be subtracted would depend strongly on their 
properties, and in particular how fast and by how much the local 
slope of the power law describing it changes. If the foregrounds are 
well-fit by a single power law over the relevant frequencies, then 
there is a good ch ance to detect the global transition signal, see 
IShaver et al.Ul999h for more detailed discussion. 

The brightness temperature behaviour we find suggests that 
the best frequency range to aim for could be ~ 120 — 150 MHz, 
since there the signal decrement is largest. A slightly higher over- 
all normalization of the primordial density fluctuation power spec- 
trum, as suggested by combining all the available data sets would 
shift the best frequency range to slightly lower frequencies, but 
ones which are still above the FM range (u > 110 MHz). Only 
a fairly high normalization, as ~ 0.9 and no power spectrum tilt 
(as given by WMAP1) would shift the global step into the FM range 
(90 MHz < v < 110 MHz). However, such a high normalization 
is largely excluded by the best current fits to the cosmological pa- 
rameters. 

We also note that albeit at first sight it appears that modifying 
the background cosmology is degenerate with assuming different 
source efficiencies, in fact the two options are not equivalent. A 
different background cosmology changes the timing of reioniza- 
tion to earlier or later, i.e. is simply a shift in redshift due to the 
delay or acceleration of st ructure formation, as previously pointed 
out in lAlvarez et al. 1 120061) . On the other hand, changing the source 
efficiencies has a somewhat different effect - either extending or 
shortening the duration of reionization - for given power spectrum 
parameters the sources form at the same time and same numbers, 
the different efficiencies just mean that their effect is stronger or 
weaker, the reby extending or sho rtening reionization, as we have 
discussed in Mellema et al. ( 2006b). While this is a somewhat sub- 
tle distinction, it does yield different shapes of dTb(z), so the two 
cases can in principle be distinguished and are thus not equivalent. 

An important reionization signature to look for is the one due 
to individual, rare, bright features. The magnitude of such features 
depends on the scale observed, ranging from ~ 0.1 K for high 
resolution (0.25', 30 kHz), to few tens of mK for beams of size 
a few arcminutes. The peak magnitude is only weakly dependent 
on redshift, with the peak value within a factor of two of the typ- 
ical values. The redshift/frequency at which the peak is reached 
is beamsize-dependent, moving to higher frequencies for larger 
beams. Taking also into account the much larger fluxes correspond- 
ing to larger beams, this again clearly argues for utilizing relatively 
large beams (5'-10' or more) and aiming at the high frequencies 
(~ 130 — 160 MHz) for such observations. The statistics of such 
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Figure 8. kSZ temperature fluctuation maps from simulations: r'250 with large-scale velocities (top left), f250C with large-scale velocities (top right), f250 
(bottom left), and f'250C (bottom right). (Images produced using the Ifrit visualization package of N. Gnedin). 



rare, bright peaks are also favourable at late times, we find that at 
~ 3' — 5' scales there are up to an order of magnitude more such 
peaks than a Gaussian statistics would predict. 



4 KSZ EFFECT FROM PATCHY REIONIZATION 

Next, we turn our attention to the second main reioniza- 
tion observable, the small-scale secondary CMB anisotropics 
gener ated through the kinetic S u nyaev-Zel'dovich (kSZ) ef- 
fect JZeldovich & Sunvaevl 1 19691: ISunyaev & Zeldovicb] Il98fj. 
Ostriker & Vishniac 1986; Vishniac 1987; Jaffe & Kamionkowski 
19981: iGruzinoy & Hul 1 1998b iHul boool: iGnedin & Jaffd l200ll ; 
Springel et ail l200ll: iMa&Frvj 120021: ISantos et alj | 2003| : 
Zhang etal.1 |2004|; Izahn et al.l 120051: iMcOuinn et al.l bood : 
Salvaterra et al.l 12005b liliev et all l2007bh . The kSZ effect results 
from Thomson scattering of the CMB photons onto electrons 
moving with a bulk velocity v. Along a line-of-sight (LOS) defined 
by a unit vector n the kSZ temperature anisotropics are given by: 

AT 



Tomb 



drje 



v ' an e OT — 



(6) 



where r\ = J dt' /a(t') is the conformal time, a is the scale factor, 
<jt = 6.65 x 10~ 25 cm -2 is the Thomson scattering cross-section, 
and r cs is the corresponding optical depth. Since our simulation 



volume is not large enough to follow the complete photon light 
path from high redshift to the end of reionization, we have to com- 
bine several simulations volumes. In order to avoid artificial effects 
from the box repetition (in particular having the same structures 
repeating over and over along the photon path) we randomly shift 
the simulation box along both directions perpendicular to the ray 
direction and also alternate the x-, y- and z-directions of the box. 
We have presented our detailed methodology along with our pre- 
dictions for the kSZ reioniza tion signal based on the WMAP1 cos- 
mology in Ilievetal. (2007b). Some of the power from the largest- 
scale velocity field perturbations is missing from the simulation 
data, due to our finite box size. We compensate for that missing 
p ower by adding it in a statistical way, again as discussed in detail 
in llliev et all J2007bh . 



4.1 The kSZ signal 

The resulting kSZ maps are shown in Figure[8] The full scale of the 
maps corresponds to our full simulation box size (100 h~ 1 Mpc, 
or ~ 50') and the pixel resolution corresponds to the simulation 
cell size (100 h~7 2 03 Mpc ^ 0.5h _1 Mpc, or ~ 0.25')- At a 
few arcminute scale there are a number of fluctuations larger than 
5 i-iK with both positive and negative sign. The f250 and f250C runs 
yield largely similar level of fluctuations, slightly larger ones in the 
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Figure 9. PDF distribution of <57ksz/TcMB (solid) vs. Gaussian distribu- 
tion with the same mean and width (dotted) with (top panels) and without 
(bottom panels) correction for the large-scale velocity power missing from 
the computational box for our simulations, as labelled. 



latter case. The typical scale of the kSZ temperature fluctuations in 
run f250 is also a bit larger than the scale for f250C, reflecting the 
larger, on average, sizes of the ionized regions in the former case 
compared to the latter. Introducing the correction for the missing 
large-scale velocity power increases the fluctuation range by about 
50% and introduces some large-scale coherent motions. 

In Figure [9] we show the PDF of these kSZ maps. The full 
range of the temperature fluctuations at pixel level is ±20 /jK 
(±13 /iK) with (without) the large-scale velocity corrections. The 
rms fluctuations of AT/Tcmb for run f250 are 4.835 x 10~ 7 
(8.968 x 10~ 7 ) without (with) large-scale velocities, while for run 
f250C the numbers are 5.094x 10~ 7 (1.020 x 10~ 6 ) without (with) 
large-scale velocities. Both the range and the rms of the kSZ fluctu- 
ations are lower than the corresponding quantities we derived from 
the WMAP1 reionization scenarios by factor of ~ 2, thus the PDF 
distributions are correspondingly less wide. In both WMAP3 cases 
there are some mild departures from Gaussianity at the bright end. 
Adding the correction for the large-scale velocities yields wider (by 
factor of ~ 2) and, interestingly, much more Gaussian PDF distri- 
butions. 

Since the kSZ effect is a product of electron density and ve- 
locity, a naive expectation is that the kSZ angular power spectrum 
would scale as <r§ for the density and a 2 for the velocity, <rf in to- 
tal. The actual ratio we find is erf' 94 for reionization scenario f250C 
and erf' 83 for f250. However since 8h 1 Mpc is well above the 
scales relevant to reionization, thus in the presence of tilt a more 
appropriate scaling would be a bandpower at dwarf galaxy scales, 
o"o.i> introduced in §[2] rather than at galaxy cluster scales, erf. We 
find the kSZ power scales as one power less, oc ctq' 74 for f250C 
and oc cr§'f 2 for f250. This scaling also holds true for the mini- 
halo bandpower, Oo.oi- Thus it becomes insensitive to the exact 
length scale once we are in the relevant part of the power spec- 
trum. As we mentioned in §[JJ the values of ero.oi and co.i scaled 
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Figure 10. Sky power spectra of STkgz/TcMB fluctuations resulting from 
our simulations: f250 (blue), and f250C (red). Solid (dashed) lines show the 
results with (without) the correction for the large-scale velocity power miss- 
ing from the box. compared to the after-reionization kSZ signals (assuming 
overlap at z ov = 8): linear Ostriker-Vishniac effect, labelled 'OV (long- 
dashed, pink), and a fully-nonlinea r model ma tched to high-resolution 
hydrodynamic simulations of Zhang et al. (2004), labelled 'Zhang et al' 
(short-dashed, dark green). We also show the total (patchy reionization + 
Zhang et al. post-reionization) signals (dot-dashed; top lines) and the pri- 
mary CMB anisotropy signal (dotted, brown). 



to include the growth inhibition between high redshift and redshift 
zero, {a/D) » Q^ 23 , 0.74 and 0.76 for the WMAP3 and WMAP1 
cases, provide a good indication when the hierarchy spanning the 
collapse of first the minihaloes then the dwarfs developed. 

The output from radio telescopes is sky power spectra, shown 
in Figure [TO] The kSZ signal from patchy reionization dominates 
the primary anisotropy at small scales, for I > 4000. The magni- 
tude of the signal is {£(£+l)C e /2Tv] ~ 10~ 13 , or AT ~ 1/xK.The 
presence of sub-grid gas clumping boosts the power by ~ 50% for 
I > 4000, but has little effect on the power spectrum shape. Adding 
the correction for the missing velocity power at large-scales boosts 
the signal power by additional factor of 2-3 in the interesting range 
of scales (t > 4000). At larger scales the boost is larger, by up 
to an order of magnitude, but at those scales the kSZ temperature 
anisotropy is strongly dominated by the primary CMB. Compared 
to the predicted post-reionization anisotropy signals, the kSZ from 
patchy reionization is larger than the linear Ostriker-Vishniac sig- 
nal (OV) for t < 20, 000 and slightly lower than, b ut similar to 
the ful l nonlinear post-reionization effect prediction by Zha ng et"al] 
(2004) (which we rescaled to the current WMAP3 cosmology us- 
ing the scaling Cg, oc g|). The tota l , p atchy reionization and post- 
reionization (based on Zhang et all d2004l) ) signals are also shown 
at the top of the figure. The total power spectrum retains the peak 
at i ~ 3000 — 4000 imprinted by the patchy reionization com- 
ponent, although the signal decrease at small scales is not as pro- 
nounced, since the decrease is partially compensated for by the 
continuing weak rise of the post-reionization component of the sig- 
nal. The two reionization scenarios (red and blue top curves) would 
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be very difficult to distinguish based solely on the total signal, but 
this might be possible to do if the post-reionization component is 
known sufficiently well and is subtracted to a good precision (see 
also the next section). A slightly higher power spectrum normal- 
ization, as suggested by recent cosmological parameter estimates, 
will moderately increase the signal as indicated by the above scal- 
ings. However, the strongly-peaked shape of the signal will persist, 
since it is reflecting the characteristic scales of the reionization pro- 
cess. Higher source efficiencies would shift the peak some towards 
larger scales, reflecting the larger typical size of the bubbles (due to 
the combination of the higher source efficiencies and the stronger 
source clustering at higher redshift) in this case. 

4.2 Observability: kSZ from patchy reionization 

The Atacama Cosmology T e lescop e jFowleJ |2004 
iKosowskv & the ACT Collaboration 2006), currently entering 
into operation, will observe at three frequency channels, at 147, 
215 and 279 GHz, targeting clear atmospheric windows, with 
bandwidths 23, 23 and 32 GHz and at resolutions 1.7', 1.3' and 
0.9', respectively. The target sensitivities in the three channels 
are 300, 500 and 700 fiK s -1 ' 2 , respectively, with a final aim of 
~ 2 (i K per pixel over a large are a of the sky (~ 200 — 400 deg 2 ). 
The South Pole Telescope (SPT) jRuhl & e t al. 2006) will be ob- 
serving in 5 frequency channels, 95, 150, 219, 274 and 345 GHz, 
with similar bandwidths and resolution to ACT. Its sensitivity 
might be even better, reaching ~ 10 fiK over 1 deg 2 in an hour. 
Thus, in terms of both resolution and sensitivity either bolometer 
is well-set to detect the patchy reionization kSZ signal. 
The thermal noise of the detectors are given by 



Ni = (s6) exp 



£(£ + l)b 2 
8 In 2 



(7) 



assuming white noise with rms s and a Gaussian beam with FWHM 
of b. The error bar corresponding to a bin Al is then given by 

2 



{ACt 



-{Ct+Nt,tatY 



(8) 



(21 + l)A£f sky 

where Ci is our patchy reionization signal, f s ^ y is the sky cov- 
erage fraction of the survey and Ne,tot = Ni + Ce )Pr i mary + 
Ce, P ost-reion is the average statistical noise for that bin. In the last 
expression we added the primary and post-reionization signals to 
the noise, since for the purposes of patchy power spectrum mea- 
surement we assume that the primary CMB fluctuations are well 
normalized, and can be statistically subtracted, contributing to the 
statistical noise. Sim ilarly, the post-re ionization kSZ signal was 
forecast robustly bv lZhang etail d2004l) and would be subtracted 
from the power spectra, but contributes to the statistical errors. 
In Figure QT] we show our predicted kSZ signal for both of our 
WMAP3 reionization models along with the ACT (left) and SPT 
(right) expected sensitivities, for s — 3.3 fj, K, b = 1.1', 100 deg 2 
area (A CT; / sky = 100/41253 = 0.0024: iHuffenberger & Seijakl 
J2005I) ) and s = 11.5 [iK, b = 1', 4000 deg 2 area (SPT; / aky = 
4000/41253 = 0.097). This assumes perfect subtraction of all 
other foregrounds. Results show that the reionization signal should 
be observable with both ACT and SPT and in principle could even 
distinguish different reionization scenarios. A number of difficult 
problems remain, however. 

The detected signal would be a mixture of thermal Sunyaev- 
Zel'dovich (tSZ) from galaxy clusters, gravitational lensing- 
induced anisotropies, Galactic dust and extragalactic point sources 
(e.g. dust in high redshift galaxies), in addition to the kSZ patchy 



reionization and post-reionization components. Separating these 
signals from each other presents significant challenges. The tSZ 
signal, which tends to dominate at these frequencies could be sep- 
arated through its characteristic spectral shape, and in particular 
using the fact that the signal goes to zero at 217 GHz in the Earth 
frame. Detecting the galaxy clusters with tSZ could allow also their 
kSZ contribution to be evaluated and subtracted, through detailed 
modelling of the clusters based on their tSZ data. The bright point 
sources could be subtracted based on complementary observations 
with e.g. Atacama Large Millimetre Array (ALMA). The lensing 
contribution is spectrally the same as the kSZ signal (and the same 
as the primary CMB anisotropies), but is statistically-different from 
the kSZ, which sh ould allow their separation, at least in principle 
dRiquelme & Spergell 20061) . 

The most difficult problem is to separate the patchy reioniza- 
tion and the post-reionization kSZ signals since both their spec- 
tra and their statistics are the same. Such separation is required, 
however if we want to extract any reionization information from 
the detected kSZ signal. This could be done for example by suf- 
ficiently detailed modelling of the post-reionization signal and its 
properties. The linear effect, also called Ostriker-Vishniac (OV) ef- 
fect, can now be calculated with a reasonable precision, but it sig- 
nificantly underestimates the actual post-reionization signal. The 
full nonlinear kSZ post-reionization effect is still difficult to derive 
from simulations due to insufficient dynamic range. Current mod- 
els and simulations roughly agree, but only within a factor of 2 at 
best, which would not allow for precise enough subtraction. An- 
other option is to use the characteristic, fairly sharp peak of the 
patchy signal at £ of a few thousand, which is in contrast to the 
much broader peak of the post-reionization signal. 



5 SUMMARY AND CONCLUSIONS 

We presented detailed predictions of the signatures of inhomoge- 
neous reionization at redshifted 21 -cm line of hydrogen and kSZ- 
induced CMB small-scale anisotropies. Our results are based on 
the largest-scale radiative transfer simulations to date, utilizing a 
background cosmology given by the WMAP 3-year data. We dis- 
cussed the observability of these signals in view of the expected pa- 
rameters and sensitivities of current and upcoming 21-cm and kSZ 
experiments. We suggested some observational strategies based on 
our results. In particular, the best approach for detecting the red- 
shifted 21-cm observations is to utilize relatively large beam sizes 
(a few arcminutes or more) and bandwidths (hundreds of kHz), 
which would result in large gains in flux, while retaining most of 
the signal. Additionally, it is better to concentrate on the high fre- 
quencies, above 120 MHz, since the 21-cm fluctuations, the corre- 
sponding fluxes and instrument sensitivities peak there, while the 
foregrounds are noticeably lower than they are at lower frequen- 
cies. 

While these basic features of the reionization signals remain 
valid for any reionization scenario in which most ionizing radiation 
is produced by stars in galaxies, the detailed figures are dependent 
on the assumptions made about the reionization parameters. These 
parameters are still highly uncertain and include the ionizing source 
efficiencies (how many ionizing photons are emitted per unit time) 
and gas dumpiness at small scales (how many recombinations oc- 
cur). It is difficult to estimate these parameters since the simulations 
do not yet have the enormous dynamic range required and since we 
do not understand sufficiently well the processes of star formation, 
galaxy formation and escape of photons from galaxies. Instead, the 
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Figure 11. Observability of the kSZ: the sky power spectrum of the reionization signal (black, solid; simulations f250 and f'250C) with the forecast error bars 
for ACT (left) and SPT (right) vs. I. The primary CMB anisotropy (dotted) and the post-reionization kSZ signal (dashed) are also shown and are added to the 
noise error bars for the reionization signal. Cosmic variance from tSZ would increase the error bars, but here we are assuming that the tSZ component has 
been completely separated, by virtue of its characteristic spectral shape. 



approach we have taken is to make simple, but reasonable assump- 
tions for these parameters and vary them within the range allowed 
by the current observational constraints. As new and more detailed 
observations become available over time, these will impose much 
more stringent constraints, which, combined with further detailed 
simulations will be able to tell us more about the properties of 
galaxies and stars at high redshifts. 

Another important caveat is that the results presented in this 
work are based on reionization simulations which do not resolve the 
smallest atomically-cooling halos, with masses from ~ 10 8 to ~ 
2 x 10 9 Mq and the even smaller moleculaiiy-cooling minihaloes. 
Smaller-box, higher-resolution radiative transfer simulatio ns which 
included all atomically-cooling halos jlliev et alj|2007al ) showed 
that the presence of low-mass sources results in self-regulation of 
the reionization process, whereby r cs is boosted, while the large- 
scale structure of reionization and the epoch of overlap are largely 
unaffected. This is a consequence of the strong suppression of these 
low-mass sources due to Jeans-mass filtering in the ionized regions. 
We expect that this self-regulation would not affect our current re- 
sults significantly since the reionization signals discussed in this 
work are dominated by the large bubbles. Utilizing smaller com- 
putational boxes in order to resolve the low-mass sources makes 
the problem more manageable, but results would underestimate the 
large-scale power of the ionization fluctuations and be subject to a 
large cosmic variance. Resolving all halos of mass ~ 1O 8 M0 or 
larger in 100 h~ x Mpc box would require ~ 10 11 particles, with 
the additional complication that on such small scales gasdynam- 
ical effects also become important and thus the complete treat- 
ment would require a fully self-consistent N-body, gasdynamics 
and radiative transfer. While still quite difficult, such simulations 
are now becoming possible with the available algorithms and com- 
puter hardware. Future higher-resolution calculations with more 
detailed microphysics would allow us to evaluate more stringently 
the effects of low-mass sources and small-scale structure on the 
reionization observables. 
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